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Abstract 

The importance of properly taking into account the factorization scheme dependence of 
parton distribution functions is emphasized. A serious error in the usual handling of this 
! topic is pointed out and the correct procedure for transforming parton distribution functions 

^ I from one factorisation scheme to another recalled. It is shown that the conventional MS and 
■ DIS definitions thereof are ill-defined due to the lack of distinction between the factorisation 

o . 

, of the strong coupling constant as- A novel definition of parton distribution functions is 

^ ! suggested and its role in the construction of consistent next-to-leading order event generators 

^ ' briefly outlined 

iu: 

> 
X 

d ' During recent years significant progress in the determination of parton distribution functions 
(p.d.f.) in the nucleon has been achieved, basically as a result of new data [|I], 0], combined 
with more sophisticated and reliable theoretical analyses P, ^ . In the CTEQ Collaboration 
^ a number of theorists, phenomenologists and experimentalists have combined their efforts 
in order to deal properly with all experimental and theoretical subtleties of quantitative QCD 
analysis of vast amount of data from various experiments and processes. In ^ p.d.f. are 
determined with high accuracy, unheard of just a few years ago. In such circumstances a 
careful reanalysis of various theoretical uncertainties is clearly needed. Although most of such 
uncertainties are discussed in sufficient detail in review papers like 0, 0], there is one which 
has not been so far satisfactorily covered in either these or in any other paper I am aware 
of. It concerns the factorization scheme (FS) dependence of finite order QCD predictions in 
processes involving hadrons in the initial state. The treatment of this ambiguity presented 
in 0, 1^, ^ is incomplete and moreover contains an error in the very central point of the 
factorization mechanism. As the FS dependence of p.d.f. has so far obtained much less 
attention than it probably deserves, I discuss in this note several of its aspects, drawing on 
analogy with the much more publicised case of the renormalization scheme (RS) ambiguity 
of the running coupling constant (couplant) a^. I think that much of the confusion and 
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misunderstanding that surrrounds this topic stems from the failure to distinguish these 
related but in principle separate uncertainties. 

The organization of this paper is as follows. In Section 2 the notation is introduced and 
some basic facts about the RS dependence of the couplant are recalled. The crucial point, 
i.e. the dependence of p.d.f. on the choice of the FS is discussed in Section 3, followed in 
Section 4 by a few critical remarks on the currently used p.d.f.. In Section 5 the merits of 
the so-called "zero" FS are discussed, and in particular it is shown how it can be used for 
the construction of consistent next-to-leading order (NLO) event generators. The results are 
summarized and conclusions drawn in Section 6. 



2 Remarks on renormalization scheme dependence 

Before coming to the ambiguity in the definition of p.d.f. let me recall a few basic facts 
about RS dependence of the renormalized couplant a(/i) = g^{fi)/4:Tr^. In massles QCD (to 
which I restrict my attention) it obeys the equation 

= -foa^ (l + ciaifi) + C2a^{fi) + ■■■), (1) 
dm fi ^ ' 

where the coefficients 6, Ci, are fixed by the number of quark flavours, while all the higher ones 
are essentially free, defining the so called renormalization convention (RC), RC={cj,i > 2}. 

In the simplest case this couplant enters the perturbation expansion of a physical quantity 
i?, depending on a single external momemtum Q, in the form 

R{Q) = a{n) (l + ri(/i/Q,ci)a(/i) + r2(/i/(5, Ci, C2)a^(/i) H ) . (2) 

Although not written out explicitly, also the couplant a(/x) (when (|lD is considered to the k-th 
order) depends on all Cj, i < k — 2. Moreover, both the couplant and the coefficients depend 
also on the specification which of the infinite number of solutions to (|I]) we have in mind. 
Each of these solutions can be labelled, for instance, by the familiar A parameter. Combining 
this last information with that on q defines what is usually called the renormalization scheme. 
Only if this RS is fixed does the specification of the scale fi (together with q, i > 2) uniquely 
determine both the couplant and the coefficients r^. Although I prefer the terminology 
advocated in P], where the the term "RS" is reserved for a unique specification of both 
a and (for detailed discussion of this point, see |^), I adopt in the following the more 
conventional notation in order to stay in close contact with ^ |^ . In this notation /i is set 
equal to some "natural" scale in the problem and the variation of the RS is parametrized by 
means of the corresponding Ars and the coefficients Ci,i > 2. Considering now (0,0) to the 
NLO, we are left with only one degree of freedom, corresponding to the variation of Ars. To 
make the following considerations as transparent as possible, let me set ci = in the rest of 
this note. [] The internal consistency of perturbation theory implies, using (|l|), the following 

^No essential conclusion obtained in the following does, however, depend on this purely technical 
simplification. 
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relation between /x, a and ri: 



riWQ,RS) = Mn(-^)-p=--p ^ a = ^— (3) 
\Ars/ a ri + p 



p = 61n(-^) -ri(/i = Q,RS) (4) 



where p 

n = /)ln I 

.Ars- 
is the renormahzation group (RG) invariant [§, which contains the Q-dependence of R{Q). 
As the variations of fi and the RS are actually two sides of the same coin, it is redundant 
to vary both of them. Without the los of generality we can fix p (for instance by setting it 
equal to some external momentum, like p = Q) in and elsewhere and vary the RS only. 
On the other hand, as the "natural" scale in the problem is usually not so unambiguously 
defined after all, many authors vary both the scale /i and the RS. Although unnecessary 
complication, this procedure is certainly legal. One has, however, to keep in mind that 
without the specification of the RS the choice of /i doesn't fix either the couplant a(/i) or the 
coefficients r^. In different RS's the same /i implies different a(/i/AR,s) and and thus the 
choice of the RS is as important as that of /i. With this in mind let me continue to label the 
RS by means of Ars, but keep p still as a free parameter. Eq. (|]) furthermore suggests that 
instead of Ars, the value of ri can equally well serve the purpose of labelling the various RS. 
Substituting into (Q) and truncating it to the NLO we get 

as an explicit function of ri and p. The obvious consequence of the nontrivial dependence 
of on the RS is that it would be a profound mistake to "transform" the couplant from 
one RS into another by means of equating the NLO (in fact any finite order) approximation 
to (0) in two different RS, say RS*^^\ RS'^^\ i.e. by solving the equation 

a(RS('))(l + ri(RS(^))a(RS('^)) =a(RS('))(l + ri(RS('V(RS('))) (6) 

Assuming (^ in RS^""^^ and expressing ri(RS^^'') in terms of a'^^^ = a(RS^^''), a*^^) = a(RS''^^) 
we get from (H) 

which yields the correct relation (^ between ri(RS*-^'') and a*^^-* only for the trivial case 
In other words imposing the relation (P) leads to inconsistency for any nontrivial 
RG transformation! In the particular case when RS*^^-* is defined by means of the effective 
charges approach of corresponding to rf^^ = 0, (|) implies 




:"7./rj (8) 

while the correct relation reads a^^^ = 
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3 Parton distributions in general factorization scheme 



In processes involving hadrons in the initial state, there is, beside the mentioned RS depen- 
dence of the couplant, another kind of ambiguity, concerning the definition of p.d.f. beyond 
the leading order [ll]]. Again, to simplify the discussion as much as possible, I restrict the dis- 
cussion to the nonsinglet (NS) quark distribution functions as revealed in the lepton-nucleon 
deep inelastic scattering (DIS). Contrary to the RS dependence of the couplant [T2|, |13 



this latter ambiguity has so far received much less attention [|T4|. This is somewhat surpris- 



ing, taking into account that its consequences may actually be even more important than 
those discussed in Section 2. 

In DIS the factorization theorem |]1T[ implies that the generic NS structure function 



F {x,Q ) can be written as a convolution (I drop the label "NS" in the following) 

F{x, Q') = f %iy, M)C f ^, -, a(/i)) (9) 
Jx y \M y ) 

of the perturbatively uncalculable quark distribution function M), defined at the fac- 
torization scale M, and obeying the evolution equation 



dlnM Jx y 



+ ■ 



(10) 



and the hard scattering cross-section C {Q / M , z , a{fj,)) admitting perturbation expansion in 
powers of the couplant at the hard scattering scale /i, generally different from M 

C{Q/M, z, a{fi)) = 6{1 -z)+ a{fi)C^^\Q/M, z) + ■ ■ ■ (11) 

While P(°) (z) is unique, both P^^^ in (|T0]) and C*-^-* in (p^ ) are ambiguous, but internal 
consistency of the factorization procedure links the variation of the NLO hard scattering 
cross-section C^^^ with that of the NLO branching function P^^^ 



C^^\Q/M, z, RS) = P(°)(z) ln(g/M) + P^^\z)/h + k{z, RS) (12) 

where the FS-invariant function n{z, RS) still, however, depends on the RS of the couplant 
a(M). Similar consistency conditions do exist at each order of perturbation theory. In the 
rest of this note I stay within the NLO approximation. The dependence of C^^^ on the RS 
of the couplant appears as a consequence of the fact, that the r.h.s. of (|l^) is given as an 
expansion in powers of this couplant. I shall return to this point in the next section. For 
fixed RS of the couplant a(M), P^^\ or via (ID, C^^\ defines at the NLO the FS of p.d.f. : 
FS:={P^^(2;)}. In the above relations the factorization scale M is again, as in the case of /x, 
kept as a free parameter. 

The first point I want to emphasize is that the renormalization of the couplant, to the 
NLO order fully described by variations of Ars in both the couplant a and the coefficients 

^If considered to all orders of a{p), C{Q/M, z, a{p)) doesn't actually depend on fi ||T^. Contrary to the 
M-dependence of q{N, M) and C^^^ {N, M), which is the basic feature of the factorization theorem and holds 
to all orders, the dependence of C{Q/M,z) on ^ is merely a consequence of truncating expansion ( pl| ) to a 
finite order. 



4 



rfc, is independent of the FS of p.d.f. as specified by the NLO branching function P^^^ in the 
evolution equation (0)! Any combination of the RS={Ars} and FS={P(^''} represents in 
principle equally legal definition of p.d.f.. Moreover, as P^^-* is a function of z, it represents 
in fact an infinite number of of degrees of freedom, and thus its variations can be expected 
to be at least as important as that of the factorization scale M. Also this point will be 
elucidated in the next section. 

Secondly, let me recall the obvious fact that although the physical observable F{x, Q"^) 
is, when are taken to all orders, independent of the parameters describing the renor- 

malization and factorization schemes, any finite order approximation to these expansions 
inevitably leads to nontrivial dependence of F{x, Q"^) on M, RS and FS. I mention it here 
as it is related to the basic question I want to address in this section and which concerns the 
way the p.d.f. transform when the FS={P''^^} is varied. As all the considerations are much 
more transparent in terms of conventional moments, defined, for a generic function f{x), as 

f{N)= [\''-'fix)dx (13) 
Jo 

let me rewrite ([9|-|T2|) in terms of them, explicitly writing out the dependence on both 
the RS of the couplant a{M) and the FS of the p.d.f. {dN,d^^\ k,{N) are the moments 
of P(°) (z) , P(i) (z) , k{z) respectively) : 

rfg(Ar,M,RS,FS) . ^ , .m^,, 



d\nM 
This is easily solved: 



g(iV, M, RS,FS)77v; 7^ = c/jva(M, RS) + d)^^(FS)a'(M, RS). (14) 



g(Ar,M,RS,FS) = A^(a(M,RS))-^^/''exp -a(M, RS)ci^^(FS)/6 , (15) 



7(1)/ 



where the constants An, introduced in [11], are independent of M as well as d^^. For 
moments of the structure function (^ we get, excplicitly writing out the dependence of the 
NLO approximation to F{N,Q'^) on M and FS={rfS^^},: 

P(A^,Q^M,RS,(iS^^) = g(iV,M, RS,rfS^) (l + a{i^,RS)C^^\Q/M, N^RS^d'ff)) , (16) 

with the following consistency condition, implied by (p!2|), 

C^'\Q/M, N, RS, 4^) =dN\n^ + ^ + k{N, RS). (17) 

For each moment A^, the expression ( [T6| ) is a function of M which, however, still depends 
on two futher parameters, one specifying the RS of the couplant and the other (d^^) the FS 
of p.d.f.. If we now want to transform g(A^, M, RS, FS) from one FS={(ij^''} into another, 

specified by FS={(i^''}, we again cannot do so by imposing the relation (the RS-dependence 
is suppressed in the rest of this section) 

q{N, M, 4^) (1 + a(M)C(^) {Q/M, N, 4^)) = q{N, M, d^) (l + a{M)C^'^ {Q/M, N, d^^ 
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As emphasized above, F{N,Q'^, M,d^^^) does nontrivially depend on M,d^j^\ as well as the 
RS and thus postulating equation like (|l^) would violate this basic feature of finite order 
approximations. Formally this is clear from inserting (|15|) into (|18D and solving the resulting 

equation. Similarly to (|^), we get the following relation between C(^) = C'^^\Q / M , N ^ S}^) 



and d!}^: 

4T = 4._jL_,Ji±^W£!il (19) 
" " a(M) \l + a(M)Cml 

which reduces to the correct one, as given in (|T^, only for the trivial case = d%^\ For 
any other case the equation (|19D is incompatible with the the consistency condition ([I^ ) 
and thereby wrong. I discuss this point in detail as in many papers, including [0, §|, 01; the 
equation (converted into moments and restricted to the NS channel) used to transform the 
quark distribution function between the so-called DIS and MS "schemes" (more on them in 
the next section) 

gDis(iV, M) = gsjs(iV, M) (l + a{M)&^{N, M)) (20) 
is (C^Vs = by definition) precisely of the incorrect form (|ll 



The only theoretically consistent way of transforming g(A^, M, FS) from one FS={(ij^^} 
into another is given explicitly in (|15D with, as emphasized, the constants held fixed. In 
I have discussed the whole procedure, based on the use of Jacobi polynomials [jl5|, |T6[, in 



X— space. Although currently other, superior, methods of solving the evolution equations are 
available [H, the fact that Jacobi polynomials are constructed from conventional moments 



(15) for which we know how the FS transformations operate, makes them invaluable in this 
kind of considerations. 

Finally a remark. The constants Ajq represent the most natural way of parametrizing 
the uncalculable nonperturbative properties of the nucleon. They are not related to any 
particular "initial" Mq, nor to any FS={(i^'*}, but determine the asymptotic behaviour of 
q{N, M, RS, d^^^) as M oo, which IS unique. 



4 Remarks on current phenomenology 

In the preceding Section I have discussed the central question of the FS dependence of the 
p.d.f.. Let me now turn to the current phenomenology of DIS, related to this subject. 

The first remark concerns the meaning of the words "MS" and "DIS" , when used in the 
connection with the p.d.f. at the NLO. The MS p.d.f. are defined in as those "which 
appear in the equation such as (^ with hard scattering part C^^^ calculated with the MS 
subtraction prescription." Although correct, this definition is obviously incomplete, as it 
specifies merely the RS of the couplant but tell us nothing about the FS={P^^^(^)} to be 
used in (^O])! Recall that the term ln47r — 7e defining the MS "subtraction scheme", is 
an artifact of extending the definition of the couplant into 4-e dimensions. As any RS of 
the couplant may be combined with any FS of p.d.f., there is an infinite set of MS-like 
p.d.f., sharing the same definition of the couplant, and therefore the right to be called "MS", 
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but arbitrarily differing in P^^\ The one usually understood in the literature under the 
label "MS" corresponds to the FS used in |18|, which should be called "MS" as it sets 
(within the dimensional regularization), the finite parts of the counterterms, renormalizing 
the appropriate composite operators, to zero. Altough the name for a FS is basically a 
matter of semantics, the current use of the term "MS" is misleading as it fails to specify the 
FS of the p.d.f. used. 

Ambiguity of a different kind is associated with the use of the label "DIS" . This scheme 
is supposed to be defined by the condition: 

cSsiQ = M, z, RS,FS) = PSk-2> RS) = -bK{z, RS) (21) 

which, however, is again not unique, due to the fact that the FS invariant k{z, RS) still 
depends on the RS of the couplant. In fact it is the the combination 

e{z) = Kiz, RS) + (z) ln(Q/ARs) (22) 



which is independent of M, FS={P^^\z)} as well as the RS of the couplant |T^. Conse- 
quently in a given RS of the couplant a(M, RS) we find 



For example we find 



pSU^, RS) = -be{z) + hP^'\z) In (23) 



PSk^, MS) = P« (^, MS) + 6PW(Z) In ^ (24) 



In general, there is again an infinite set of "DIS"-like p.d.f., which share the property C^^-* = 
0, but differ in the NLO branching function P'^^\z) and the couplant, thereby leading to 
different numerical predictions when inserted into (^). P.d.f. bearing the name "DIS", like 
those of |T3, or some of [ffl, tacitly assume MS as the RS of the couplant. This, however. 



is not a must and thus for an unambiguous specification of the "DIS" -like FS the RS of the 
couplant should always be specified. 

The second remark concerns the practical aspect of exploiting the vast freedom in the 
definition of p.d.f. at the NLO. As already mentioned, little phenomenological attention has 
so far been payed to the FS dependence of p.d.f.. This may be due in part to the failure to 
appreciate the independence of these two renormalization procedures. It is also true that to 
apply, for instance, the idea of "optimization" P to the FC dependence of p.d.f. in x-space is 



technically much more involved. In [|T4| I have, however, argued that at least the FS defined 
by setting P^^^ = (at the NLO; at higher orders it would generalize by setting all higher 
order AP branching functions to zero) should seriously be considered. In this "zero" FS the 
full NLO correction is put into the hard scattering cross-section C^^\ thereby representing 
in some sense the opposite of the DIS FS, which sets C^^^{z) = 0. Moreover, it turns out 



1^ that when the moments of p.d.f. are considered, this FS is very close to that obtained 
via the Principle of Minimum Sensitivity of |^ (for ci = they even coincide). Although 
this results doesn't automatically imply the same close relation for the p.d.f. themselves, it 
seems reasonable to add this FS to the list of those used in phenomenological applications. 
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5 "Zero" FS and NLO event generators 

There is, in fact, another reason why this FS could be of considerable interest. Recall, that all 
currently used event generators, like HERWIG, PYTHIA, JETSET, LEPTO etc., are based 
on essentially leading-log parton showers. Although they are sometimes combined with NLO 
hard scattering cross-sections, the overall description remains only LO. The simple picture 
of LO parton showers becomes much more complicated when one attempts to generalize 
them to the NLO. To get a consistent NLO description of any hard scattering process in the 
"zero" FS we, however, need merely the LO parton showers as P^^\z) = in this FS! As the 
example of moments of structure functions shows, this choice may be quite reasonable and 
should definitely be tried. Although simple at first glance, one has to be careful in taking 



for the NLO cross-section that corresponding to this "zero" FS. From ([12D we easily find its 
form: 

Ci'jM/M, z) ^ PW(^) HQ/M) + Kiz, RS) (25) 

Using the results of p[ on P^^) and C^^\ k{z, RS) can straightforwardly be evaluated in the 
MS RS. Transformation to any other RS is then trivial. 

Let me stress that this procedure is not equivalent to the so-called matching of parton 
showers to fixed order matrix elements, as recently implemented in LEPTO event generator 
| |2CI| |. There, the exact 0{as) matrix element is matched to the parton shower at some 
particular value of incoming parton virtuality tm (see Fig.l) in the sense that below tm only 
parton showers are used while above tm the matrix element takes fully over. In the case of 
"zero" FS, the situation is different and the NLO cross- section da^^^{Q/M, z,t)/dt which, 
when integrated over t, yields (^Sj) contributes at any virtuality, even below that given by the 
factorization scale t = M^, as only the pole term 1/t plus some finite part is factorized into 
the parton p.d.f.. As a result, the NLO hard scattering cross-section da^^^ {Q /M, z,t) / dt 
becomes a discontinuous function of t at t = M^, this discontinuity being cancelled by a 
similar discontinuity of the parton shower contribution, which is restricted by definition to 
the domain t < M^! Work on practical implementation of this idea is in progress. 

The situation is sketched in Fig.2a, where the full da^^'-^{Q/M, z,t)/dt, corresponding 
to the sum of the diagram in Fig.l and the one with the gluon radiated off the outgoing 
quark, is plotted as a function of r = — t for fixed and z. The full result is a sum of three 
terms: the pole term of the form A{z)/t, the residue of which is proportional to P^'^\z), 
the T independent constant B{z,Q) and C{z,Q)t, linearly rising with r. Notice that for 
r ^ 0, the finite part of dcr^^'^(z, Q, t), i.e. the sum of the last two terms, becomes negative. 
The result of factorization, i.e. the separation of the full NLO cross-section, containing all 
infrared and parallel singularities, into a part included in the quark distribution function and 
the remaining, finite, hard scattering cross-section, is represented in Fig. 2b by the dashed 
and dotted curves, discontinuous at the factorization scale r = M^. For r > the 
hard scattering cross-section da'^^'-^{z,Q,t) coincides with the full result, but for r < 
its definition is ambiguous as it depends on how much of the finite part will accompany 
the sigular pole pole term into the definition of the quark distribution function. Fig. 2b 
corresponds to the case that only the pole term A/t is subtracted. Recall that, for instance, 
the term P^°\z) \n{Q/M) appearing in (|l^) is essentially the integral of the pole term A/t 
from to the upper kinematically allowed value of r, proportional to Q^. 
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6 Summary and conclusions 



In this note I have discussed several aspects of factorization scheme dependence of parton 
distribution functions. I have emphasized potential importance of proper treatment of this 
ambiguity for the reliabhty of theoretical analyses of ever better data. Special attention has 
been payed to the correct transformation of p.d.f. between different factorization schemes 
and the ambiguities in the meaning of some of the currently most popular definition of p.d.f. 
have been brought to light. Finally the so-called "zero" FS has been proposed and shown 
to be potentialy useful in the construction of NLO event generators. 
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Figure captions 

Fig.l: Feynman diagram describing the process e~-|-q— >e~+q+g with parallel singularity in 
the ^-channel. 

Fig. 2a: A typical shape of da{Q, z,t)/dt as a function of r = — t for fixed x,Q^. In this 
example the pole term (dashed curve) is added to a linearly rising finite part (dotted line) 
to give the full NLO contribution (solid curve). 

Fig. 2b: Separation of the full NLO contribution (sohd curve) into the part absorbed in the 
quark distribution function (dashed curve) and the finite NLO hard scattering cross-section 
(dotted curve) . The discontinuities of the last two curves at r = cancel in the sum. 
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This figure "fig2-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9305240vl 



This figure "fig3-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9305240vl 



